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Abstract 

We perform the renormalization group analysis on the dynamical symmetry 
breaking under strong external magnetic field, studied recently by Gusynin, 
Miransky and Shovkovy. We find that any attractive four-Fermi interac- 
tion becomes strong in the low energy, thus leading to dynamical symmetry 
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breaking. When the four-Fermi interaction is absent, the (3- function for the 
electromagnetic coupling vanishes in the leading order in 1/N. By solving the 
Schwinger-Dyson equation for the fermion propagator, we show that in 1/N 
expansion, for any electromagnetic coupling, dynamical symmetry breaking 
occurs due to the presence of Landau energy gap by the external magnetic 
field. 
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Recently dynamical symmetry breaking in the presence of the external magnetic field 
has attracted much interest. Gusynin, Miransky and Shovkovy |I| have shown that constant 
magnetic field acts as a strong catalyst of dymamical symmetry breaking, which leads to 
fermion mass generation in 2+1 as well as in 3+1 dimension. They also have shown that 
similar phenomena happen in 3+1 QED without extra four fermion coupling [0. In a recent 
article || Dunne and Hall discussed the problem in the context of a non-uniform magnetic 
field in 2+1 dimesnsion. 

In this paper we perform the renormalization group analysis and show that the 4-fermi 
interaction becomes strong in the infrared and leads to instability that causes dynamical 
symmetry breaking. For QED under strong external magnetic fields, we do not have a 
RG interpretation for dynamical mass generation, because the electromagnetic coupling is 
not running. Thus, one has to rely on nonperturbative methods to see dynamical mass 
generation. We solve the Schwinger-Dyson equations for the fermion propagator in the 
1/N expansion and find nontrivial and energetically preferable solutions. By comparing 
the solutions with operator product expansion for the fermion two-point function, we argue 
that for any small electromagnetic coupling fermions get dynamical mass. This result is 
consistent with the analysis of Bethe-Salpeter equation for fermion and anti-fermion bound 
state done by Gusynin et al. [0. 

We first consider a free fermion under constant magnetic field in z direction, B = Bz. 
The Lagrangian is given by 

C = VipV (1) 

where If = 7 M (<9 M — ieA^) and \I/ is a massless fermion. We choose the symmetric gauge 
A = (—^y, ^x, 0) and solve the eigenvalue equation, 

[a ■ (p - eA)]V = EV (2) 

to get the spectrum and the basis of this system. The eigenvalues are indexed by collective 
index A = (a, j3, n, k z ) and given by 
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E A = ay/k* + 2 \eB\n (3) 

where a — ± denotes the sign of the energy, (3 = ±| is the spin component along the 
magnetic field, and the quantum number n is given by 

2n = 2n r + 1 + |m L | - sign(e.B)(m L + 2(3). (4) 

n is a nonnegative integer that labels the Landau level. Here n r is the number of nodes of 
radial eigenfuction, uil is the angular momentum of the eigen function. The eigen function 
is 



U A = NA ^k z z el m L <P r \m L \ L \rn L \ ^ eB ^ ^ 



(5) 



where N A is the nomalization, L r n {x) is the associated Laguerre polynomial, and u a ^ = 
Xa <8> Vp 1S the eigenvector of a 3 <E> 03 where two o" 3 correspond to the energy and the spin. 
It is a spinor in the Lorentz frame in which the magnetic field is specified. This basis forms 
an orthonormal system. The eigenfunction expansion 

*(x,t)= fj2Mt)U A (x), (6) 

J A 

yields 

So = fY,Mt)\idt-EA)Mt) (7) 

J A,t 

Finally we take the Fourier transform in t to get 

So = f d 4 x^i 0* = /X>1(^ - E A )^ A (u) (8) 



where 

A,w a,f3,n r ,m L 

We now determine the scaling dimesions of various modes. Notice that n = 0, the lowest 
mode, has a different scaling property from the rest of the modes. Under the scaling, 



k z — y sk z , uj — > su, with s < 1 (10) 

we require that the kinetic term be invariant. To meet this, if)o(k z ,u) must have scaling 
dimension —3/2, while ip n> o(k z , uS) must have —1, since the Landau gap should not scale. 
Namely, 

il> {k s ,u) -> s- 3/2 ij (u,k z ) (11) 

x/j n (k z ,u) — > s~ x ij} n {u),kz) for n > 0. (12) 

We see that the n = mode is more relevant than the other modes. This is a manifestation 
of the dominance of the lowest Landau level (LLL) or decoupling of the fermions at higher 
Landau levels. We now determine the scaling dimension of the interaction. Again, we split 
the interaction term Si nt into two parts; one contains only zero modes, the other contains 
at least one n^O mode; Si nt = S® nt + S' int , where 

s int = J n^dw < <y(^A^)5(^o; < )V'S^o^^iflor»V'o,oV'o,r> (13) 



with 



C AB cd = J dxU{(x)U B (x)uUx)U D (x) (14) 

where each of A,B,C,D contains the LLL index n = 0. The dimensionality of Sf nt is 
4(2 + (—3/2) x 2 x 2) - 1 — 1 = 0. Since the zeroth mode is the most relevant term, all 
the rest terms have positive scaling dimension therefore they are irrelevant interactions. 
Therefore the transverse directions are completely decoupled from the low energy dynamics. 
The interactions are at most marginal at tree level. This is the proof of the dimensional 
reduction which was also noticed in |1[ from the structure of the Schwinger propagator. So 
to settle the issue whether there is any relevant interaction that leads phase transition, we 
have to go to the one loop (3 function. Unless there is a special reason, the one loop (3 
function does not vanish in general. The crucial factor then is its sign. 

Consider NJL model in (3 + 1)D with external electromagnetic field, y4® xt (x), described 
by Lagrangian density, 



G r 



c = m jp^ + — 



(^) 2 + (^7 5 ^)' 



(15) 



where p = 7 /i (<9 At — ieA^) and \I/ is a massless A^-column 4-component spinor. Let A be 
some scale of the system under study. We choose it to be less than the Landau level gap 
^/jaB[. we decompose the field \l/ into slow and fast components 



qr = -qj + -qj 



f» 



(16) 



where 



E 



- - / e ikzZ ^(x,y,k z ), 

k z \<sA 



*/= e ik * z V{x,y,h z 

sA<\k z \<A 



(17) 



Since the transverse momentum should not scale we do not split it into slow and fast. We 
integrate the fast modes out and see how much change is induced in the action of slow part. 
For example out of (G/2) 2 [( 1 l' s + 'i> / )(\I> s + \I> / )] 4 we are looking for the coefficient of (\i> s \I> s ) 2 . 
We only get 2! x 2! times the fish diagram. Since the calculation is straightforward we only 
mention the structure of the propagator. 

As derived by Schwinger |3[], the fermion propagator in a constant external magnetic 
field is 



S(x,y) = S(x-y) exp 



(x- y yA°*(x + y) 



(18) 



We take the magnetic field in z direction. The Fourier transform of 5* is given as 
with 



n=0 



kl - kl -2 \eB\n 



(19) 



D n (eB,k) = (k° 1 °-k 3 1 3 ) (l-* 7 Vsign(e5))L n (^ 



- 1 + ry 7 sign(e£?) ) L 



■>n-\ 



2k\ 
\eB\ 



+ Ak ± ■ 7j.Li_! 



\eB\ 

2k\ 
]eW\ 



(20) 
(21) 



where L" are the generalized Laguerre polynomials and L n = l? n , = 0. k± is momentum 
perpendicular to the external magnetic field. As we have shown above, only LLL fluctuation 



contributes to the low energy dynamics so we can use the propagator that is projected to 
LLL instead of the Schwinger's full propagator. 

§(k) = l e~Al\^\h P _ + o ( , (22) 



kf ' " \\eB\, 

where k u = k — k± and P_ = 1 — i7 1 7 2 sign(ei?) is the projection operator to project out 
fermions of spin parallel to the external magnetic field if e > or antiparallel if e < 0. The 
one loop correction to the four-Fermi interaction is 

z6 G = (f ) 2 iV f sA<|fc|<A Tr [s(k)S(k) + S(k^S(k^] (23) 
= 2G 2 iV/ sA<|fc|||<A e- fc i/l^l- fc i/l^l™P 2 (24) 
= ^\eB\G 2 h\\ (25) 

If we set g = \eB\G with dimensionless coupling g, then 5g = \eB\SG. Hence the above 
equation can be written as 

The dynamical mass generated in the infrared can be determined by integrating the 
renormalization group equation Eq. fl26|). We find 



^dyn ~ V \eB\e ^d^i) (27) 

The magnitude of dynamical mass is determined by the strength of the four-Fermi coupling 
at scale yfaBj. This result agrees with the vacuum energy calculation of Gusynin et al. 

Now, let us consider the case when four- Fermi interaction is absent, namely, pure QED, 
in which case the symmetry of the system is enlarged; 

SU(N) — > SU(N) L x SU(N) R (28) 

In general, the electromagnetic coupling is weak in the IR region. Therefore nonperturbative 
effects like dynamical symmetry breaking do not occur unless QED is in the strong phase 
0]. But, under a strong external magnetic field, since the excitations of fermions occur only 
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in the magnetic field direction, fermions may interact strongly. To see this one has to rely 
on another nonperturbative method. Here we use 1/N expansion, keeping e 2 N = a finite 
when N — > oo. Then, for the photon propagator, we have to sum all bubble diagrams. The 
one loop vacuum polarization at low energy is 

U af3 {p) = -{-ief -N^Ty \s{k) 1 a S{p + k)^] (29) 
= +a J k e-^-^k^Tr [|p_ 7 <* j^P-^] (30) 
= ^\eB\e-vl/^ B \(gf- P yP/pi), (31) 

where g^f = g a/3 for a, (3 — 0, 3 otherwise zero. The photon propagator in 1/N expansion is 
then in Landau gauge 

where 

n{p) = pL\eB\ e -^^. (33) 

We see that the singularity of the photon propagator in the parallel direction is shifted from 
p = to nonzero p. Photon develops mass along the parallel direction m 2 = ^ \eB\, while 
it remains massless along the perpendicular direction. The Coulomb potential between two 
static fermions is 

/oo p t4 -mr 

J X ° J j2^ W ° 0{p)etP ' X ~ (34) 

Since the vacuum polarization tensor U a ^(p) is finite, the electromagnetic coupling is not 
renormalized at the leading order in 1/N expansion. 

The simplest order parameter for symmetry breaking is the condensate of fermion bilin- 
ear; 

(M) = lim^ \x - y\ lm (V(x)y(y)) (35) 
= (V^o) , (36) 



where 7 m is the anomalous mass dimension of ty^?(x) and in the second line we used the 
fact that fermions at higher Landau levels decouple. One can extract the value for the order 
parameter from the operator product expansion of the fermion propagator |J. Solving the 
Schwinger-Dyson equations, we determine the fermion propagator. For energy less than the 
Landau gap, y|aB[, the full fermion propagator will be of following form; 

g W=fe -*l/| eB| _l_ p,. (37) 

Note that the fermion self energy Z and £ are independent of k±, since the on-shell condition 
has to be 

k* + m 2 phy = 0, (38) 

where m p ^ y is the physical mass. Schwinger-Dyson equation for the fermion propagator in 
the leading order in 1/N is 

[ ( z W - 1) A - e«] p_ = | S^'^ z^t-nk/ ^- (39) 

At leading order, we can take Z(p u ) = 1. Then the (chirally non-invariant) fermion self 
energy is, after Wick rotation and taking trace over Dirac matrices, 

Nj (JJtt)* 6 kl + E(& n ) 2 (p n - A;,,) 2 + k\ ' [W) 

After integrating over k± and isolating the logarithmic divergence at (p u — k u ) 2 = 0, we 
expand the r.h.s. of Eq. ( f40"D in powers of (p„ — k u ) 2 / \eB\ and get 

a f E(jfe n ) r 1 



ln ^' -0.063 + 

in \eB\ 



(41) 



N J h kl + £(fc„) 2 

where the ellipsis denotes terms of order of (p u — k u ) 2 / \eB\ or higher which are negligible. 
(From now on, we let p = \p u \ and k = \k u \.) After differentiating Eq. (|4l| ) with respect to 
p, we integrate over the direction of k n to get 

d v a f p „k E(k) , , 



Multiplying by p and differentiating once again, we obtain 



pYl 1 + T! + r- 



pY, 



0. 



p 2 + Y? 

where the prime denotes differentiation with respect to p and we have defined 



(43) 



a 



2n 2 N' 



(44) 



An infrared boundary condition is given by Eq. fl4"2|); 



limpS'(p) = 0. 



(45) 



One can easily find the soulutions to the second-order, nonlinear differential equation ( |4"3"| ) 
in two asymptotic regions. For p <C S(p), 



Y,{;>) ~ //><• -f A In ( - ) . 



(46) 



But, because of the IR boundary condition ([E|), A = and S(p) ~ mc for p <C S(p). For 
p > E(p), 



S(p) ~ m 



(47) 



j? / \p 

where e = i^fr. The parameter /z, with mass dimension corresponds to the renormalization 
point. By using same analysis done by Cohen and Georgi H, one can show that there are 
solutions which connect two asymptotic solutions. 
We let 



X{t)=Y(p)/p, t = \n{p/ti). 



(48) 



Then Eq. (H3) becomes 



X + 2X 



dX 



V(X), 



(49) 



where 



V(X) = -[X 2 + rHl + X 2 )} 



(50) 
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This is the equation for a particle of unit mass moving in a potential V under a friction 
proportional to veleocity. The IR boundary condition implies that as t —>■ — oo the particle 
is sitting at 

X(t) - ^e~*. (51) 

Since, for large X, V ~ |X 2 , the particle moves toward X = 0, critically damped by the 
friction. But, near X = 0, the particle is underdamped. Therefore the particle sitting at 
X(— oo) will eventually get to X = in a finite time and then oscillate around X = 0, 
corresponding to the asymptotic solutions (|^). We see therefore there must exist solutions 
which interpolate the IR and the UV solutions. Futhermore, since the the vacuum energy 
at the leading order in 1/N is 



V(E) 



eB\ f d 2 k 



4tt J (2tt) 



-2 \n(k 2 + S 2 ) + 2 



E 



2 



(52) 



such solutions have lower vacuum energy than the trivial solution, S = 0. 

The operator product expansion for the fermions in the lowest Landau level has the 
following form: 



lim <^o(p)> =^+^ + --- (53) 

— p- 1 — >oo " y 

= ^+p(f) < + ^(f)" < + -, (54) 

where p 2 = w 2 — k 2 z . (Note that two operator m and ip^ip^x) have same mass dimension 
but opposite anomalous mass dimension.) But, since the nontrivial solutions to the gap 
equation (|43|) have oscillatory behavior in deep UV region (due to imaginary anomalous mass 
dimension), one cannot tell from the operator product expansion of the fermion propagator 
which is the renormalized mass and which is the fermion condensate. But, as a signal for the 
dynamical symmetry breaking, we see the renormalized mass and the fermion condensate 
coalesce for any value of electromagnetic coupling a. Therefore, as a general rule, we have 
dynamical symmetry breaking, (^Z 1 ^) ^ in QED in a strong magnetic field. 

In conclusion, we analyzed dynamical symmetry breaking under external magnetic field in 
the Wilsonian renormalization picture. We found that any attractive four- Fermi interaction 
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is enhanced in the infrared, leading to dynamical mass generation for fermions. When the 
four- Fermi interaction is absent, the /3-function for the electromagnetic coupling is zero in 
the leading order in 1/N. By solving the Schwinger-Dyson equation, we found that any 
small electromagnetic coupling leads to dynamical symmetry breaking, just because at low 
energy (E < y/\eB\) fermions at high Landau levels decouple and fermions at lowest Landau 
level are (l+l)-dimensional and photons are propagating in (3+l)-dimensions. It will be 
of great interest to extend this result to the non-homogeneous magnetic field context. As 
a prelimanary step to this direction, one can study the small deviation from the constant 
magnetic field. In this case one can easily show that this deviation contains at least one 
power of momentun and it is irrelevant interaction term, therefore does not chage the ground 
state qualitatively. This background deviation, of course, is different from the quantum 
fluctuation of the photon field whose effect is shown to be relevant in this paper indirectly. 
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